Calorons, instantons and constituent monopoles in SU(3) lattice gauge theory 
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We analyze the zero-modes of the Dirac operator in quenched SU(3) gauge configurations at non- 
zero temperature and compare periodic and anti-periodic temporal boundary conditions for the 
fcrmions. It is demonstrated that for the different boundary conditions often the modes are localized 
at different space-time points and have different sizes. Our observations are consistent with patterns 
expected for Kraan - van Baal solutions of the classical Yang-Mills equations. These solutions consist 
of constituent monopoles and the zero-modes are localized on different constituents for different 
boundary conditions. Our findings indicate that the excitations of the QCD vacuum are more 
structured than simple instanton-like lumps. 

PACS numbers: 11.15.Ha 



Understanding the relevant excitations of the QCD 
vacuum is a major scientific goal and lattice calculations 
can contribute valuable non-perturbative information. A 
prominent candidate for such relevant excitations are in- 
stantons which are particularly successful in providing 
a mechanism for chiral symmetry breaking [1]. Lattice 
calculations have helped to corroborate these ideas. Re- 
cently several studies of the zero-modes and near zero- 
modes of the Dirac operator [2-4] which provide a pow- 
erful and natural filter retaining the infrared fluctuations 
of the underlying gauge field were presented. It was es- 
tablished that the vacuum has a lumpy structure and the 
lumps are locally chiral [3] and (anti) self-dual [4]. 

However, it is generally believed, that instantons are 
not responsible for confinement. On the other hand all 
lattice calculations show that the deconfinement transi- 
tion and the restoration of chiral symmetry coincide at a 
single critical temperature which suggests that there is an 
underlying mechanism linking the two phenomena. In re- 
cent years Kraan and van Baal [5] have constructed new 
solutions for the classical SU(N) Yang Mills equations 
with compactified time, i.e. at finite temperature (KvB 
solutions) . The KvB solutions generalize the caloron so- 
lution [6] and a lump of topological charge 1 is built of N 
constituent monopoles, a feature which might provide a 
link to the confinement problem. Note that for KvB solu- 
tions the monopole lumps appear in multiples of N. A sin- 
gle lump is not an independent object and only N lumps 
together have topological charge ±1. Strong evidence for 
SU(2) KvB solutions in cooled lattice gauge configura- 
tions were given for twisted [7] and periodic boundary 
conditions [8]. The zero- mode for the KvB solution [9] 
shows characteristic features which should allow to use 
the eigenvectors of the Dirac operator as a filter retain- 
ing infrared fluctuations and search for traces of KvB 
solutions also in thermalized configurations. 

Before presenting our lattice calculations we briefly dis- 
cuss a few properties of SU(3) KvB solutions and their 
zero-mode. KvB solutions [5] are characterized by the 
Polyakov loop at spatial infinity which in a suitable gauge 



can be expressed as Poo = exp(i27rdiag(^i, /i2, M3)) with 
Mi + + Ma = and \i\ < [i 2 < /13 < M4 = 1 + Ml- The 
action density can be written down in a relatively simple 
form [5] and depends in addition to the phases ^i,/i<2,M3 
also on three 3- vectors y\, j/2, 2/3- As remarked the KvB 
solutions can be seen to be superpositions of constituent 
monopoles and the 3 corresponding lumps are in space 
located at y%, ij2, 2/3. The m-th constituent monopole can 
be assigned a mass 8ir 2 (fi m+ i — fi m ) (both the masses 
and the monopole positions y\ are in units of the inverse 
temperature) . A general choice for the parameters of the 
KvB solution will give rise to three lumps at different 
positions with different width. 

The zero-mode i/j in the background of a KvB solution 
has been computed [9] for arbitrary temporal boundary 
conditions ij)(t + l/T,x) — exp(i2nz)ip(t, x) with z — 1/2 
corresponding to anti-periodic boundary conditions (b.c.) 
and z — giving the periodic case. The remarkable 
feature of the zero-mode is that it traces only one of 
the monopoles. In particular it is located on the m-th 
monopole when z € {/i m ,/i m+ i}. At least for well sepa- 
rated monopoles the width of the lump seen by the zero- 
mode depends on the size of the underlying gluonic lump 
and also the quantity a = min[z — /i m ,/i m+ i — z]. The 
lump is narrow for large a and broad for small a. This 
is an important feature in the deconfined phase where 
due to the Z3 symmetry of the action the Polyakov loop 
can have three different phases 0,±27r/3. At sufficiently 
high temperatures configurations with real Polyakov loop 
have \ii ~ 0, % = 1, 2, 3. The anti-periodic zero-mode has 
z = 1/2, is located on monopole 3 and has a ~ 1/2 
and thus is quite localized. The periodic zero-mode with 
2 = can be located on any of the monopoles. It has 
a ~ and is very broad. For configurations with phase 
27r/3 of the Polyakov loop (the case with phase — 27r/3 
is similar) one has (p,±, fi2, M3) ~ ( — 2/3, 1/3, 1/3). The 
anti-periodic mode has a ~ 1/6 while the periodic mode 
has ~ 1/3 and both of them are localized on monopole 
1. This implies that for complex Polyakov loop the roles 
are reversed and here the periodic mode is more local- 
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statistics 


a [fm] 


T [MeV] 


T/T c 


6 x 20 3 


8.20 


70 


0.115 


287 


0.96 


6 x 20 3 


8.45 


64 + 25 


0.094 


350 


1.17 



TABLE I. Parameters of our calculations. 



ized. However, the difference in localization is not as 
pronounced since the values a ~ 1/6 and 1/3 differ less 
than a ~ 1/2 and obtained for real Polyakov loop. 

To summarize, if KvB solutions are present in thermal- 
ized SU(3) gauge configurations then one should find in- 
stances where the zero-modes with anti-periodic b.c. and 
periodic b.c. are located at different positions and the 
lumps they see can differ in size. In the dcconfincd phase 
the differences in size should show a particular pattern 
when comparing configurations with real and complex 
Polyakov loop. These are simple stable signatures which 
distinguish KvB solutions from instantons or calorons. 

In our numerical analysis we use quenched SU(3) gauge 
configurations generated with the Liischer-Weisz action 
[10] on 6 x 20 3 lattices. We adjust the inverse gauge 
coupling (3 such that one ensemble is below the critical 
temperature and the other one in the deconfined high 
temperature phase. The lattice spacing a and the val- 
ues for the temperature T and T c for the Liischer-Weisz 
action listed in Table I were determined in [11]. 

We calculate the eigenvalues and eigenvectors of the 
chirally improved Dirac operator [12] using the implicitly 
restarted Arnoldi method [13]. This is done for both, pe- 
riodic and anti-periodic temporal b.c. for the fermions 
(the spatial b.c. are always periodic). The chirally im- 
proved Dirac operator is not exactly chiral but a good 
approximation of an exactly chiral lattice Dirac opera- 
tor. The topological zero-modes of the continuum mani- 
fest themselves as small real modes. It can be shown [14] 
that although the eigenvalues are not exactly zero they 
give rise to the correct topological susceptibility through 
the index theorem. Although the small real eigenvalues 
do not vanish exactly we will refer to the correspond- 
ing modes as zero-modes in the following. From a larger 
set of configurations we select those configurations which 
have a single zero-mode for both b.c, i.e. are in the sec- 
tor with topological charge Q = ±1. The restriction to 
the |Q| = 1 sector is necessary to avoid mixing between 
multiple zero-modes that occurs in higher sectors. For 
our configurations with \Q\ = 1 we find that for both, 
periodic and anti-periodic boundary conditions the zero- 
mode is always dominated by a single lump. The position 
of the lump, however, can be different for the two choices 
of the boundary conditions. 

The ensemble in the high temperature phase has a 
non-vanishing Polyakov loop. We will use this ensem- 
ble to study the influence of the expectation value of 
the Polyakov loop on the zero-modes. We divide the 
deconfining ensemble into sub-ensembles with complex 
respectively real Polyakov loop containing 64 (25) con- 




FIG. 1. Scalar density of the zero-mode for anti-periodic 
temporal b.c. (l.h.s.) and periodic temporal b.c. (r.h.s.). We 
show y-z slices for a configuration below T c . 



figurations. We will refer to these subensembles as the 
complex respectively the real sector. 

In order to analyze the localization of the zero modes 
we form the scalar density p{x) = ^ c d \ip(x) c ,d\ 2 from 
the zero-modes ip{x) c ^d by summing over the color and 
Dirac indices c, d. The resulting density p is gauge in- 
variant. In Fig. 1 we show y-z slices of the scalar density 
and compare the anti-periodic case (l.h.s.) to the peri- 
odic case (r.h.s.). The plots are for a configuration in the 
confining phase. 

Fig. 1 clearly demonstrates that the zero-modes with 
periodic and anti-periodic b.c. are located at different 
space-time points. For anti-periodic b.c. p has its maxi- 
mum at (t, x, y, z) — (2, 9, 13, 13) while for periodic b.c. it 
is located at (1,9,6,7). The figure shows the y-z slices 
through the maxima, i.e. the t = 2, x — 9 slice for anti- 
periodic b.c, respectively the t = l,x = 9 slice for the 
periodic case. The maxima are separated by 9.3 lattice 
spacings which corresponds to 1.1 fm. It is also obvious 
from the figure that the two lumps have different height 
and width. This is not what is expected for the zero- 
mode of a classical caloron but agrees well with the be- 
havior predicted for the KvB solutions. When scanning 
through our ensemble we have detected many configura- 
tions showing such behavior. 

Having found instances of configurations resembling 
KvB solutions we now turn to a more quantitative anal- 
ysis. One interesting observable was already addressed, 
the distance between the maximum of pa and the maxi- 
mum of pp. The subscript A refers to p computed from 
the zero mode with anti-periodic b.c, while the subscript 
P is used for p computed from the periodic zero-mode. 

A related observable is the overlap of the support of 
the two lumps. The support is defined as the function 
9(x) which has 6{x) = 1 for the lattice points x carrying 
the lump in p and 9{x) = for all other points. An 
open question is how to identify the lump. For a narrow 
lump the support is small while a broad lump has a much 
larger support. 

We will base our definition of the support on the so- 
called inverse participation ratio / = V ^ x p(x) 2 , where 
V is the total number of points in our lattice. Note that 
we use normalized eigenvectors such that ^2 X p(x) = 1. It 
is easy to see that for an entirely localized state (p(x) = 1 
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FIG. 2. Histograms for the distance of the maxima (l.h.s.) 
and the relative overlap of the lumps for different boundary 
conditions (r.h.s.). The data are for the ensemble below T c . 



for only one lattice point and else) the inverse partici- 
pation ratio is I — V , while for a completely delocalized 
state (p(x) = 1/V for all lattice points) 1=1. The 
inverse participation ratio is a convenient measure for 
the localization, assuming large values for localized states 
and small values for spread out states. 

Using / wc now define a number N = V/ (32/) and 
use this to determine the support by setting 9(x) = 1 for 
those ./V lattice points x which carry the N largest values 
of p and 9(x) =0 for all other points. The normalization 
of N was chosen such that if the lump was a 4-D Gaussian 
the support would consist of the 4-volumc carrying the 
points inside the half-width of the Gaussian. Note that 
N is the size of the support, i.e. N — J2 X 0(x). 

With our definition of the support 9 for both the 
anti-periodic (9 a) an d periodic (Op) lumps we can de- 
fine the overlap function C of the two lumps as C (x) = 
9p(x)9a(x). From this we compute the relative overlap 
R = 2/(N A + N P ) J2 X C(x). The relative overlap R is a 
real number ranging from 0.0 to 1.0. If the lump in p A 
and the lump in pp sit on top of each other (at least the 
core up to about the half-width) and have the same size 
then R = 1.0. If the two lumps are entirely separated 
one has R = 0.0. Values in between give the ratio of the 
common volume of the two supports to the average size 
of the two supports. Note that even if the two lumps sit 
on top of each other, but have different sizes Na, Np one 
has R < 1.0. 

In Fig. 2 we show histograms for the distance between 
the maxima in p A and pp (l.h.s.) and histograms for 
the relative overlap (r.h.s). The distances were rounded 
to the nearest integer while the relative overlap R was 
rounded to the nearest multiple of 0.1. The data shown 
are for the ensemble below T c . When looking at the 
l.h.s. plot one finds that for many configurations the max- 
ima of the two lumps sit on top of each other. However, 
the distribution shows also that for 50% of the configu- 
rations the two lumps are separated by at least 6 lattice 
spacings which corresponds to 0.7 fm, almost twice the 
typical instanton size. This indicates that for a large 
fraction of the ensemble the zero-modes with different 
b.c. see different lumps. 

Also from the histograms for the relative overlap 
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FIG. 3. Same as Fig. 2 but for the ensemble above T c . 



(r.h.s. plot in Fig. 2) one comes to the same conclusion. 
Only for a small fraction of the configurations one finds 
that the lumps seen by anti-periodic and periodic zero- 
modes coincide and for many cases they do not overlap 
at all. A comparison of the two plots in Fig. 2 shows 
that although many pairs of lumps may have the same 
position of their maximum they still differ in size such 
that their relative overlap is smaller than 1. 

In order to further study the evidence for KvB so- 
lutions wc now turn to the ensemble in the high- 
temperature phase. For these configurations the theory 
is deconfined and the Polyakov loop has a non-vanishing 
expectation value. As discussed above the phase of the 
Polyakov loop plays an important role for the properties 
of the KvB solutions and one expects significant differ- 
ences between the complex and the real sector. 

Let us begin our discussion of the deconfined phase by 
again looking at the histograms for the distance between 
the peaks of p and their overlap. Fig. 3 is equivalent to 
Fig. 2 but now shows data for the ensemble above T c . In 
addition to showing the data for the full ensemble we also 
display the results for the subensemble in the real sector. 
These data are represented by thinner, black bars. 

When inspecting the histograms for the distance be- 
tween the peaks (l.h.s. plot) we find that again a certain 
amount of configurations shows rather large distances but 
the fraction of these configurations is smaller. 12% of 
the ensemble show distances of at least 7 lattice spac- 
ings (again corresponding to ~ 0.7 fm). Note that these 
configurations are entirely dominated by the real sector. 
This is exactly as pointed out in our summary of the 
KvB zero-modes where we discussed that in the complex 
sector the two zero-modes sit on the same lump. 

Another effect is evident from the two plots in Fig. 3. 
Although more than half of the configurations have the 
center of the two lumps at exactly the same lattice point 
we do not see many configurations where the lumps have 
a large overlap. The distribution is instead peaked at 
an overlap of 0.5. The reason for this behavior is the 
difference in size between the two lumps. One of them is 
much narrower than the other one and thus their relative 
overlap is considerably smaller than 1. 

We directly study the phenomenon of the large dif- 
ference of the sizes of the two lumps by comparing the 
inverse participation ratio / of the zero-mode for anti- 



3 



200 

> 

150 
100 
50 





• real sector 

v complex sector 



•• • • 



50 100 150 200 



and a substructure resembling the constituent monopolcs 
of Kraan - van Baal solutions is evident. 
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FIG. 4. Comparison of the inverse participation ratio of 
the zero mode for periodic (Ip) and anti-periodic (Ia) b.c. 
The data are from the ensemble above T c . 



periodic and periodic b.c. (Ia,Ip)- Fig. 4 shows a scat- 
ter plot of our results with the data in the real sec- 
tor represented by filled circles and open triangles are 
used for the complex sector. In the real sector one finds 
that the anti-periodic inverse participation ratio (ranging 
from I a ~ 25 up to I a ~ 180) is much larger than its 
periodic counterpart (Ip ~ 2-5). For the complex sec- 
tor the situation is similar but the role of the boundary 
conditions is reversed. However, the effect is not as dra- 
matic as in the real sector. For not too localized modes 
(Ip up to 70) one finds that the inverse participation ra- 
tio for the periodic b.c. is about 3-times larger than for 
the anti-periodic case. The qualitative pattern for the 
different sectors of the Polyakov loop matches the pre- 
dictions for KvB solutions discussed in our summary of 
the properties of the KvB zero-modes. 

We remark that we have performed an equivalent anal- 
ysis also on a 16 4 lattice at a = 0.094 fm. Although 
this is a typical lattice for zero-temperature calculations 
we still find a substantial amount of configurations with 
zero-modes showing the characteristic properties of KvB 
solutions. This observation might account for the failure 
of identifying the profile of classical instanton zero-modes 
in quenched lattice configurations [15]. 

In this letter we have studied zero-modes of the Dirac 
operator for gauge configurations with topological charge 
±1 for temperatures slightly below T c and a second en- 
semble in the deconfined phase above T c . Analyzing the 
localization of the zero-modes we observed large differ- 
ences when comparing temporal anti-periodic and peri- 
odic b.c. for the fermions. For a substantial portion of 
the ensemble we found that the zero-mode is located at 
different positions for different b.c. and also the width of 
the lumps can differ substantially. For the ensemble in 
the deconfined phase we observed that the behavior also 
depends on the phase of the Polyakov loop. The emerging 
picture for the zero-modes is qualitatively well described 
by the properties of the zero-modes of Kraan - van Baal 
solutions. Our results indicate that the picture of topo- 
logical excitations being single lumps might be too simple 
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